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ABSTRACT 

We present results from monitoring the multi-waveband flux, linear polariza- 
tion, and parsec-scale structure of the quasar PKS 1510—089, concentrating on 
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eight major 7-ray flares that occurred during the interval 2009.0-2009.5. The 
7-ray peaks were essentially simultaneous with maxima at optical wavelengths, 
although the flux ratio of the two wavebands varied by an order of magnitude. 
The optical polarization vector rotated by 720° during a 5-day period encom- 
passing six of these flares. This culminated in a very bright, ~ 1 day, optical 
and 7-ray flare as a bright knot of emission passed through the highest-intensity, 
stationary feature (the "core") seen in 43 GHz Very Long Baseline Array im- 
ages. The knot continued to propagate down the jet at an apparent speed of 
22c and emit strongly at 7-ray energies as a months-long X-ray/radio outburst 
intensifled. We interpret these events as the result of the knot following a spiral 
path through a mainly toroidal magnetic fleld pattern in the acceleration and 
collimation zone of the jet, after which it passes through a standing shock in the 
43 GHz core and then continues downstream. In this picture, the rapid 7-ray 
flares result from scattering of infrared seed photons from a relatively slow sheath 
of the jet as well as from optical synchrotron radiation in the faster spine. The 
2006-2009.7 radio and X-ray flux variations are correlated at very high signifl- 
cance; we conclude that the X-rays are mainly from inverse Compton scattering 
of infrared seed photons by 20-40 MeV electrons. 

Subject headings: quasars: individual (PKS 1510—089) — polarization — gamma 
rays: general — radio continuum: galaxies — X-rays: galaxies 



1. Introduction 



Attempts to understand relativistic jets of blazars have been greatly advanced by the 
availability of instruments such as the Very Long Baseline Array (VLBA), the Rossi X-ray 
Timing Explorer (RXTE), and, most recently, the Fermi Gamma-ray Space Telescope, to- 
gether with more traditional telescopes. Long-term multi-waveband monitoring of a number 
of blazars wit h these facihties is now providing valuable insights into the physica l processes in 
the jets (e.g., iMarscher et al.ll2008l : IChatterjee et al.ll2008t iLarionov et al.ll2008l ). This paper 
presents a rich set of observations of the 7-ray bright quasar PKS 1510—089 [z = 0.361), 
whose jet exhibits apparent motions of er nission features that are among the fastest (as high 
as 45c) of all blazars observed thus far (jjorstad et al.ll2005[ ). We analyze motions of fea- 
tures in the parsec-scale radio jet alongside variability of the optical polarization and flux 
from radio through 7-ray frequencies. The relative timing of correlated variations probes 
the structure and physics of the innermost jet regions where the flow is accelerated and 
collimated, and where the emitting electrons are energized. 
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2. Observations and Data Analysis 



Our observations include imaging with the VLBA at 43 GHz in both total and linearly 
polarized intensity, with angular resolution near 0.10 milliarcseconds (mas), or a projected 
distance of 0.50 pc (for Ho = 71 km s~^ Mpc~^ and the concordance cosmology: ISpergel et al. 



20071). We proce s sed t he data and created images in a manner identical to that described 
bv ljorstad et all J2005h . 



We derived 0.1-200 GeV 7-ray fluxes by analyzing data from the Large Area Telescope 



(LAT) of the Fermi Gamma-ray Space Telescope with the standard software ( lAtwood et al. 



2OO9I ). Photon counts from a circular region of radius 20° centered on PKS 1510—089 were 
fit by single power-law spectral models of this source, PKS 1502+106, and PKS 1508—055, 
plus model v02 of the Galactic and extragalactic backgrounds. The mean slope of the photon 
spectrum fo r 2009 .0-2009.5 for 7-day integrations, —2.48 ± 0.05, matched that reported by 
Abdo et al.l ( 120091 ) for 2008.6-2008.83. Our final fiuxes are from 1-day integrati ons using this 



slope , with a detection criterion that the maximum-likelihood test statistic (IMattox et al. 
1996h exceed 9.0. 



We obta ined 2.4-10 keV X-ray fluxes with the RXTE PGA. We processed the data as 
described by iMarscher et al.l ( 120081 ). fitting the photon count spectrum with a single power 
law plus photoelectric absorption corresponding to a neutral hydrogen column density of 
8 X 10^° cm~^. In addition, we measured fiux densities at: 14.5 GHz with the 26 m antenna 
of the University of Michigan Radio Astronomy Observatory (see [A 
with the 14 m telescope of the Metsahovi Radio Observatory (see 



and frequencies near 230 GHz with the Submillimeter Array (see iGurwell et al 



ler et al. 1985): 37 GHz 



Terasranta et al.lll998[) : 



200 



We measured the degree of optical linear polarization P and its position angle x with: 
the 0.4 m telescope of St. Petersburg State University; the 0.7 m telescope at the Grimean 
Astrophysical Observatory; the PRISM camera on the Lowell Observatory 1.83 m Perkins 
Telescope; the Steward Observatory 2.3 and 1.54 m telescopes 0; and the 2.2 m telescope 
at Galar Alto Observatory, under the M APGAT program. The d a ta analysis pro c edure s 
for the variou s telesc opes are described in lDArcangelo et al.l ( 120071 ): iLarionov et al.l ( 120081 ): 



Jorstad et al.l ( 120101 ). We also obtained optical (R-band) fiux densities from photometric 
observations at the six telescopes hsted above, the 2.0 m Liverpool Telescope, and the 0.7 m 
Meniscus Telescope of Abastumani Astrophysical Observatory. We have added data from 



Yale University B 



^Data listing: http://james.as.arizona.edu/~psinith/Fermi 
^Data hsting: http://www.astro.yale.edu/smarts/glast 
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Observational Results and Discussion 



A number of flares are apparent in the 2008-0 9 radio to 7-ray light curv es (Fig. [H which 
includes 7-ray data from AGILE in 2008 March; iD'Ammando et al.ll2009l ). Of particular 
interest is the first half of 2009 (Fig. [2]), during which eight major 7-ray flares are apparent. 
The VLBA images (Fig. [3]) feature a bright "core," presumed stationary, from which knots of 
emission separate at apparent superluminal speeds. We detect two new knots, with apparent 
speeds of 24±2 c and 21.6 ±0.6 c. The first passed the core on JD 2454675 ±20, as the X-ray 
flux reached a sharp peak, and ~ 2 weeks after a broad maximum in the 14.5 GHz flux (see 
Fig. [T]). The second knot passed the core on JD 2454959 ± 4, essentially simultaneous with 
the extremely sharp, high-amplitude optical/7-ray flare 8 on JD 2454962 (Fig. [2]). We argue 
below that this knot was responsible for flares 3-8 and possibly 1-2. 

Flare 8 coincides with the end of a 50-day rotation of x by 720° that started near flare 3 
(see Fig. HI). This striking phenomenon can be explained either by a stochastic process or by 
a coherent magnetic field geom etry of the flaring region. In the s t ochastic interpreta tion, the 
magnetic field is turbulent (see lJoneslll988l : iMarscher et al.lll992l : lD'Arcangeldl2009l ) and the 
apparent rotation results from a random walk of the resultant polarization vector direction 
as cells with random magnetic field orien tations enter and then e xit the emission region. 
According to our simulations (described by lD'Arcangelo et al.ll2007l ). one rotation by > 720° 
similar to that of PKS 1510—089 and lasting 50 ±10 days occurs once per ~ 2, 000 days. The 
probability that such a random apparent rotation would coincide so closely with a specific 
50-day period of elevated 7-ray fiux is extremely small, ~ 0.1%. Furthermore, the model 
predicts equally probable clockwise and counterclockwise rotations of x the same object, 
while the ~ 180° rotation of x between JD 2454990 and 2455000 (see Fig. S]) proceeded at 
the same counterclockwise rate as during the final stage of the larger rotation. This implies 
that the rotation is dictated by geometry, with the accompanying minor flare occurring in 
the same location as flare 8. 



We observed a similar , but shorter, rotation of the optical polarization in BL Lac in 
2005 (IMarscher et al.ll2008l ). We apply our phenomenological model developed to explain 
that event to PKS 1510—089. A moving emission feature follows a spiral path as it propa- 
gates through the toroid al magnetic fleld of the acceleration and coUimation zone of the jet 
( Komissarov et al.ll2007l ). The spiral motion is caused either by rotation of the flow (jVlahakis 
2006[) or a helic aFstream of electron-positron pairs injected from the black hole's ergosphere 
(jWilliamsll2004j ) . The secular increase in the rate of rotation (see Fig. |4]) is due to an increas- 
ing Doppler factor as the flow accelerates until it reaches the core, since the Lorentz factor 



increases as the cross-s ectional radius of the jet (IVlahakis fc Konigll |2004J : IVlahakid 12006 



Komissarov et al.l 120071 ). [However, the rest-frame angular velocity of the emission feature 
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cann o t decrease to c onserve specific angular momentum as the jet expands (ISteffen et al. 
19951 : I Vlahakid 120061 ) . hence the spiral path followed by the centroid of the feature must 
maintain a constant radius of helical motion to avoid this.] The emission feature must cover 
most of the cross-section of the jet in order to create substantial flares and to cancel most of 
the polarization (from different orientations of the toroidal magnetic field across the feature) 
while leaving a residual of < 10%, with x rotating systematically as the feature proceeds 
down the jet. Front-to-back light-travel time delays also cancel some polarization by stretch- 
ing the feature along the spiral path as viewed by the observer. Synchrotron flares occur 
when the energization of electrons increases suddenly over some or all of the emission region; 
our model does not attempt to explain these sudden surges in particle acceleration. The 
solid curve in the bottom panel of Figure H] shows the fit to the variation of x of the model, 
with parameter values that are typical but probably not unique. The bulk Lorentz factor 
increases linearly with longitudinal distance down the jet from a value of 8 at the onset of 
the rotation to 24 at the end. The viewing angle of the jet axis is 1.4° in the model, so the 
Doppler factor increases from 15 to 38 during this time, and the final apparent speed is 21c, 
matching the observed value within the uncertainties. The infle ctions of the curve are due t o 
changing aberration, which we calculate with the equations of iKonigl fc Choudhuril ( 1l985l ). 
The scenario proposed in that paper, in which the jet twists from helical instabilities, can 
explain the rotation of x but not the observed low, randomly changing degree of polarization. 
The secular increase both in the rate of rotation of x ^i-iid in optical flux during the first half 
of the outburst results from the acceleration of the flow and consequent enhanced beaming. 



Flare 8, featuring the highest optical flux observed since 1948 ( iLiller fc Lillerl Il975l ). 
occurred as the rotation of x was ending and the new superluminal knot was passing the 
core. This event can be explained by compression of the knot by a standing conical shock 
( iMarscher et al.l 120081 ). According to our model, when the knot passed the core it was 
propagating down the jet at 0.3 pc day"^ in our frame, and had traveled 17 pc downstream 
since the rotation of x began. The transverse radius of the jet at this point i s ~ 2 x 10^^ 
cm, based on the opening angle of the jet ~ 0.2° derived by IJorstad et al.l (|2005[ ). The 



propagation rate, plus the 1-day timescale of the flare, sets the longitudinal size of the 
optical flare region at ~ 0.3 pc. We associate this size with the distance moved by an 
electron radiating at R-band before it loses too much energy to continue doing so. To 
determine the energy loss rate for all of the events, we first calculate the ratio of 7-ray 
(> 100 MeV) to synchrotron luminosity (which equals the ratio of inverse- Compton to 
synchrotron loss rate), Cgs, of the peak of flares 1-8 to be 70, 30, 40, 40, 30, 10, 40, and 
9, respectively. We estimate the luminosity of the synchrotron radiation, whose spectral 
energy distribution peaks at infrared wavelengths, as the R-band flux density multiplied by 
vr = 4.7 X 10^^ Hz, by 1.27 to correct for extinction (jSchlegel et al.lll998[ ). and then by 
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6 to convert roughly to bolometric luminosity (see iKataoka et al.l 120081 ) . We then derive 
that the magnetic field 5 ~ [(6 x 10^ s)(r/20)c/(0.3^ pc)]2/3(5/40)^/3 _ o.4^-2/3 G during 
fiare 8, where ^ = (mb + ^iphot)/^^B5 ^Cgs ~ 9 if the seed photons originate from outside the 

1/2 

emission feature, ^ ~ 1 for mainly first-order SSC, and ^ ~ ~ 3 for combined first- 
and second-order SSC. We therefore obtain B ~ 0.1-0.4 G, with the highest value for the 
first-order SSC case. The (toroidal) magnetic field at the start of the rotation of x was then 
~ 1 G. 

The 50-day rotation of x implies that a single moving emission feature was responsible for 
the entire outburst encompassing fiares 3-8. We identify this feature as the superluminal knot 
seen later in the VLBA images (Fig. [3]), whose polarization vector lies in the same direction 
as in the optical shortly after fiare 8. Our observations therefore demonstrate that the high- 
energy emission from the jet of PKS 1510—089 is quite complex, arising from different regions 
and probably by multiple emission mechanisms as a single disturbance propagates down the 
jet. Both the 7-ray and optical emission is highly variable, but not always in unison. This 
is unexpected, since when B ~ 0.1-1 G, electrons with the same energies, 7e ~ X0^'^^° '^ in 
rest-mass units, should be involved in optical synchrotron radiation and inverse Compton 
scattering of optical or IR photons to 7-ray energies. 

When the emission feature was close to the base of the jet, electrons in the jet could 
have scattered the broad emission-line region or accretion-disk photons to 7-ray energies 
as the disturbance first became optically thin to photon-photon pair production (see, e.g.. 



Ghisellini fc Tavecchid 120091 : iDermer et al.l l2009l ) . We speculate that this could correspond 



to fiares 1-2. In this case, the subsequent quiescent period would imply that there is a section 
of the jet where neither seed photons nor electrons with 7e > 10^ are abundant. 

We conclude that some or all of 7-ray fiares 1-4 and 7, with very weak optical coun- 
terparts, were caused by sudden increases in the local seed photon field at optical or IR 
wavelengths rather than by increased energization of electrons. The time scales of variabil- 
ity limit the size of each source of seed photons to a < 1 pc — not much larger than the 
cross-sectional radius of the jet, ~ 0.1 pc — while the luminosity of each would need to be 
~ 3 X 10^3(Cgs/60)(r/20)-2(5/0.4 G)2(a/0.1 pc)2erg s'^ if the source lies at the periphery of 
the jet, and higher if more remote. This luminosity is too high for any commonly occurring 
cosmic object located parsecs from the central engine, but could be obtained in a relatively 
slow she ath surrounding the u ltra-fast spine of the jet responsible for the high superluminal 



motion (IGhisellini et al.ll2005l ). Moving knots or standing shocks in the sheath could produce 
the requisite number of seed photons while being too poorly beamed to contribute substan- 
tially to the observed fiux. The relatively slow motion and gradual evolution of such features 
implies that they should persist for years, in which case another series of fiares in the near 
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future should exhibit a similar pattern of variability and appearance of a superluminal knot. 

Flares 5 and 8 included rapid optical flaring that required sudden energization of elec- 
trons to 7 ~ XO^'^''^^'^. This suggests that the flaring 7-rays could have been created by the 
SSC process, an inference that is supported by the lower value of Cgs than for the other flares. 
This ratio was, however, much greater than unit y, which implies that seco nd-order scattering 



contributed significantly to the 7-ray flux (e.g., iBloom fc Marschen 119961 ). We should there- 
fore expect a different slope of the 7-ray photon spectrum for these flares than for flares 1-4 
and 7. There are a sufficient number of photon counts to measure the slopes for flares 5, 
7, and 8: -2.36 ± 0.09, -2.73 ± 0.17, and -1.92 ± 0.21, respectively. The probability that 
these are the same is ~ 0.2%. Detailed modeling is needed to determine the slopes expected 
for the different radiation mechanisms. 

The extremely high optical flux and polarization of flare 8 imply that the magnetic 
field was higher and more ordered for this event than for other fiares. This is consistent 
with our interpretation that the 43 GHz core — where fiare 8 occurred — is a standing shock. 
The low value of Cgs implies that the magnetic field was compressed by a factor ~ 2 over 
the previous fiare. We suggest that the emission region responsible for fiare 5, when the 
polarization increased to > 20%, was a subset of the entire disturbance, so that there was less 
cancellation of the polarization. This indicates that the fiare was caused by unusually efficient 
energization of electrons over only part of the disturbance, although we have insufficient 
information to determine the cause. In contrast, fiares 6 and 7 occurred when the polarization 
was quite low, which implies that electrons were accelerated less intensely throughout the 
entire disturbance. 

The fiare in 2008 September ( JD 2454700-2454750; Fig. [I]) is of a different nature than 
flares 1-8. It took place one month after a very bright superluminal knot passed through 
th e core (which coincid ed wit h the peak of a major X- r ay flar e), consistent with the flnding 



of IJorstad et al.l ( 1200 ll ) and iLahteenmaki &: Valtaojal (120031 ) that high 7-ray states often 



follow the time when a new superluminal knot coincides with the core or a millimeter-wave 
outburst starts. The frequency- dependent delays — the emission peaked flrst at optical, 
then 7-ray (10-day lag), and then radio and X-ray (24-day lag) frequencies — contrast with 
the simultaneity of the optical and 7-ray maxima of flares 2, 5, 7, and 8. The time delays 
imply that a power law is maintained over only a limited range of electron energies, with the 
range changing as a knot separates from the core. 

Since 1996, the X-ray emission has correlated better with the 14.5 GHz variations than 



with those at higher frequencies (IMarscherl |2006| ). During the 2006-2009 period covered 



by Figure O the correlation is signiflcant at > 99% confldence. The electrons radiating 
at 14.5 GHz near the core have energies 7e ~ 20-40 in rest mass units, far too small to 
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scatter the 14 . 5 GH z photons to X-ray energies. For this reason, we conclude, as have 



Kataoka et al.l ( 120081 ). that these electrons produce the X-rays via inverse Compton scattering 



of IR seed photons from either the sheath of the jet or a source outside the jet, such as a 
dust torus (IBlazejowski et al.ll2000l ). In this case, the X-ray variability is in response to 
both the evolution of the electron energy distribution — which affects the radio emission in 
a similar way — and changes in the seed photon density as a knot propagates down the 
jet. In the observer's frame, the radiative energy loss time scale of 7e ~ 30 electrons, 
^loss ~ 3 X 10'^B~'^(~J^6~^ s, ~ 5 days when B ~ 0.4 G and 6Qs ~ 400. This is commensurate 
with the fastest X-ray and radio variability that we have observed. 



4. Conclusions 

Our comprehensive dataset has revealed some patterns that allow us to probe the loca- 
tions and causes of the high-energy emission in PKS 1510—089. The rotation of the optical 
polarization vector over 50 days in 2009, as multiple 7-ray and optical flares took place, im- 
plies that a single emission feature — the superlumin al knot later seen in ou r VLBA images — 



was responsible for these variations. As for BL Lac (iMarscher et al.ll2008[ ). we model this as 
a structure that followed a spiral path through a helical magnetic field where the jet flow ac- 
celerates. The 7-ray to synchrotron flux ratio varied greatly among the different flares. This 
requires local sources of seed photons from both within the jet and just outside, probably in 
a surrounding sheath. As a consequence, flares erupt at a variety of locations as disturbances 
pass down the jet. Future comprehensive multi- waveband monitoring plus VLBA imaging 
will determine the extent to which our findings apply to the general population of blazars. 
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Fig. 1. — Multi-waveband light curves from 2008.0 to 2009.6. Upper limits (red) are 2- 
a. R-band flux is uncorrected for reddening. Vertical arrows: times when superluminal 
knots passed the 43 GHz core; horizontal bars: uncertainties in these times. For reference, 
JD 2454700=2008 August 21, 2454900=2009 March 9. 
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Fig. 2. — 7-ray and optical light curves from 2009.0 to 2009.38. See caption to Fig. [T] 
Vertical dotted lines denote 7-ray flares discussed in text. 
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Fig. 3. — Sequence of 43 GHz VLBA images showing ejections of two superluminal knots 
with proper motions of (top) 1.1 ± 0.1, (bottom) 0.97 ± 0.06 mas yr~^. Images are con- 
volved with a circular Gaussian beam of FWHM=0.1 mas (shaded circle on bottom right). 
Calendar dates and JD— 2450000 of images are given. Contours: total intensity, with lev- 
els (top) 1, 2, 4,..., 64, 96% of peak of 1.41 Jy beam'^ and {bottom) 0.25, 0.5, 1, 2,..., 
64, 96% of peak of 3.14 Jy beam~^. White line segments: direction of linear polarization; 
color: polarized intensity relative to peak (green) of [top) 71 mJy beam~^ and (bottom) 
120 mJy be am~^. Polarization not correction for Faraday rotation, whose value is poorly 
constrained (jjorstad et al.l 120071 ). 



- 14 - 



T 



I R-Band Flux Density t 




4900 4950 5000 

JD - 2450000 



Fig. 4. — R-band flux density and degree P and electric-vector position angle x of optical 
polarization in early 2009. Filled black circles: R-band; filled triangles: V-band; open 
squares: A =500-700 nm. Multiples of 180° are added to x as needed to minimize jumps in 
consecutive values of x oi") after JD 2454990, so the values of x can be compared with the 
end of the first rotation. The curve fits the x data with the model discussed in the text. 
Highest amplitude optical fiare peaked on JD 2454962=2009 May 10. 
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Fig. 5.— Variation of X-ray, radio, and optical flux of PKS 1510-089 from 2006.0 to 2009.7. 
See caption to Fig. [TJ 



